Dynamics of the peel front and the nature of acoustic emission during peehng of an 

adhesive tape 
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We investigate the peel front dynamics and acoustic emission of an adhesive tape within the con- 
text of a recent model by including an additional dissipative energy that mimics bursts of acoustic 
signals. We find that the nature of the peeling front can vary from smooth to stuck-peeled con- 
figuration depending on the values of dissipation coefficient, inertia of the roller, mass of the tape. 
Interestingly, we find that the distribution of AE bursts shows a power law statistics with two scaling 
regimes with increasing pull velocity as observed in experiments. In these regimes, the stuck-peeled 
configuration is similar to the 'edge of peeling' reminiscent of a system driven to a critical state. 
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The process of peeling of an adhesive material from a 
substrate is a complicated phenomenon involving molec- 
ular attraction at the interface, and kinetic and dynami- 
cal effects. The kinetic nature of the process is clear from 
the fact that the peel force depends on the peel rate. The 
fracture process during peeling can be either cohesive or 
adhesive at low and high peel velocities respectively. At 
intermediate velocities, the fracture process is intermit- 
tent suggesting that the peeling process results from an 
interplay of time scales. At low peel velocities, there is 
sufficient time for viscoelastic glue to fully relax while at 
high velocities, the glue essentially behaves like a solid 
0. The intermittent behavior is observed when the vis- 
coelastic time scale is of the same order as the peel rate 
time scale. It is an everyday experience that the peel- 
ing process is always accompanied by a characteristic 
audible noise 0, 1^. However, the mechanism leading 
to the acoustic emission (AE) has remained ill under- 
stood. Moreover, the inhomogeneous deformation of the 
peel front results from the destabilization of uniformly 
advancing peel front To the best of our knowledge, 
we are not aware of any model that investigates the dy- 
namics of the peeling front and the associated acoustic 
emission. We address these two issues within the context 
of a model for the peeling of an adhesive tape. 

Experiments on peeling of an adhesive tape mounted 
on a roller show that the peel force function has two 
stable branches separated by an unstable one. The pull 
force exhibits a rich variety of behavior ranging from saw- 
tooth to irregular waveforms 2, ^] . A dynamical anal- 
ysis of the force waveforms and the AE signals reports 
chaotic dynamics at the upper end of pull velocities 0- 
However, as there are no models, no further insight into 
the origin of acoustic emission has been possible. 

A relevant model introduced in Q has been studied by 
others (Ref. ^ 5, 6J) belongs to the category of differen- 
tial algebraic equations (DAE) and are singular requiring 
an appropriate DAE algorithm provided in Ref. ^ . [For 
this reason, the results in Ref. 5] are the artifact of the 
method followed.] Recently, we modified these equations 
into a set of ordinary differential equations (ODE) by 



including the kinetic energy of the stretched tape. The 
ODE model not only supports dynamical jumps across 
the two stable branches, it displays a rich dynamics j^. 
Here, we extend this model to include spatial degrees of 
freedom to study the contact line dynamics of the peel- 
ing front. The inclusion of a local strain rate dependent 
Rayleigh dissipation functional along with the kinetic en- 
ergy of the tape forms a basis for converting the potential 
energy stored in the stretched tape into kinetic energy 
provides a mechanism for explaining qualitative experi- 
mental features on acoustic emission (AE) [kI ITiI| . 

Figure^a) shows a schematic representation of the ex- 
perimental setup. An adhesive roller tape of radius R is 
mounted on an axis passing through O and is pulled at a 
constant speed by a couple meter motor positioned at 
O' . Then, the line PQ represents the peeling front. Sev- 
eral features of the setup can be explained by considering 
the projection on to the plane of the paper {OPO'). The 
tangent to the contact point P (representing the con- 
tact line PQ) subtends an angle to the line PO' . Let 
the distance between O to O' be / and the peeled length 
of the tape PO' be L. If P subtends an angle a at O 
with the horizontal OO' , the geometry of the setup gives 
L cos 9 — —I sin a and L sin 9 — I cosa — R. As the local 
velocity v at P undergoes rapid bursts during rupture, 
we have V — v + ii — L — v + u + R cos 9 a. Let u{y) to 
be displacement with respect to the uniform 'stuck' peel 
front and let v(jj),9(y),a(jj) be defined at every point y 
along the contact line. As the tape of width b is pulled 
with a velocity V , the above equation generalizes to 
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[V - v{y) - u{y) - R cos %) a{y)]dy = 0. (1) 



However, as we are interested in the dynamics of the con- 
tact line of the softer glue material (whose elastic con- 
stant is three orders less than that of tape material) , the 
effective spring constant kg of the contact line is assumed 
to be much less than that of the tape material kt . Then, 
as the entire tape is pulled with a velocity V, the force 
along PO' equilibrates fast, we can assume that the in- 
tegrand in Eq. ^ is zero for all y. 
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We derive the equations of motion of the contact Hne 
by considering the Lagrangian L — Uk — Up, where Uk 
and U p are the kinetic and potential energies respec- 
tively. The kinetic energy is given hy Uk — \ /q^ C [0^(2/) + 

■^^j^dy + i p [«(?/)] ^(iy, where the first term rep- 
resents the rotational kinetic energy of the roller tape 
and second term arises due to the kinetic energy of the 
stretched part of the tape. Here, ^ is the moment of in- 
ertia per unit width of the roller tape and p the mass per 
unit width of the length L. The total potential energy 

(PE) Up oi the stretched ribbon can be written asUp = 

2 

dy. The total dissi- 



dy 



pation has two contributions TZ 
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dy, where f{v, V) physically represents 

the peel force function assumed to be derivable from a 
potential function = J f{v)dv (see Ref. 0). The 

second term, denoted by Rae, represents the dissipation 
arising from the rapid movement of peel front is given by 
the Rayleigh dissipative functional. This has the same 
form as the energy dissipated in the form of acoustic 
emission during abrupt motion of dislocations in plas- 
tic deformation, i.e., Eg^ oc e^(r), where e{r) is the local 
plastic strain rate \Uj^ . Hence we interpret Rgf. as the en- 
ergy dissipated in the form of AE signals. Indeed, such 
a term has been successfully used to explain several fea- 
tures of AE signals in martensites • 

We rewrite all the energy terms in a scaled form using 
a time like variable r = ujut where = kt/{h p) . Let 
fmax and Vmax the maximum value of f{v) and v on the 
left stable branch. Then, defining a length scale d — 
fmax/h, we introduce u = Xd = X^fmax/h), I = l^d, 
L = L^d and R = R^d. The peel force / can be written as 
/ = fmax4>{v^) where v" = v/vctOud and V = V/vcUJud 
are the dimensionless peel and pull velocities respectively. 
Here, Vc — Vmax/^ud is dimensionless critical velocity at 
which the unstable branch starts (in the scaled units, 
the unstable branch begins at = 1). Defining Cf = 
ifmax/ktfip/Oi h = kg9/{kta'^), 7^ = T^Uu/ {kto^), 
and y = ar, where a is a unit length variable along the 
peel front, the scaled local form of Eq. (Q) is 



X = {V' - v')vc + i?" -[^ [sin a) a (2) 



The scaled kinetic and potential energies can be written 



as Ut = 



K 2Cf Jo 
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Vf,v^ (r) 



dr+\!'J^ X{r) 



dr 



and C/|, = lj'J'^X\r)dr + !fr^ 
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dr 



dr respec- 



tively. The total dissipation in a scaled form is TZ^ = 

2 



b/a 



dX(r) 



dr 



dr. 



(w*) is 



the scaled peel force that can be obtained by using in 
Eq. (9) of Ref. [3 shown in Fig. ^h). We shall refer the 
left branch AB as the 'stuck state' and the high velocity 
branch CD as the peeled state. 

Using {a{r),d{r), X{r), X{r)) as generalized coordi- 




FIG. 1; (a) A schematic plot of the experimental setup, (b) 
Plots of (hiv") as a function of v" for = 1.45. 



nates in the Lagrange equations of motion, ^ 
0, and A ( 
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we get the equations of motion as 



dX{r) ) dX(r) ~^ 0X{r) 



X = 



R' 



CfR' 



ko 



V /L'^ sin a 
(1 + l^'/L" sin a) 
d^X <j>{v'') 



(3) 



dr"^ {\ + / L'^ sin a) 



Equations (|3I4|) are still not suitable for further analysis 
as they have to satisfy the constraint equation Eq. (|2l . In 
the spirit of mechanical systems with constraints 12], we 
obtain the equation for the acceleration variable v'^{r) by 
differentiating Eq. to be 

y-^[^X + ^(a\cosa - RH^'-^f) + sinaa)] K.(5) 
ij ij 

Eqs. (|2I3|I and © were solved by discretizing and us- 
ing an adaptive step size stiff differential equations solver 
(MATLAB package) for open boundary conditions. The 
initial conditions were chosen from the stuck state [ ie., 
AB branch of 4'{v'')\ with a small spatial inhomogeneity 
in X that satisfies Eq.© approximately. The system is 
evolved till a steady state is reached before the data is 
accumulated. We have studied the dynamics over a wide 
range of values of C/, , and 7„ keeping i?'* = 0.35, 
P = 3.5, fco = 0.1, iV = 50 and = 100. Note that Vc 
is an important parameter which however is determined 
once f{v) is given. (The values of the unsealed parame- 
ters, for example kt ~ 1000, are fixed using the data in 
Ref. f{v) used here preserves major the features 

of the experimental curve such as fmax ^ 280N/m at 
Vmax = 0.05cm/ s with a velocity jump to 16 cm/s. See 
also 0. Note however, we do not use the dynamization 
scheme used in @.) These equations exhibit rich dy- 
namics which can be classified as uniform, rugged and 
stuck-peeled nature of the contact line. In the unsealed 
variables, the results reported here correspond to chang- 
ing the tape mass (m = pb) keeping I{— ^b) constant. 

First consider the results for Cf = 0.00788 {vc — 
0.0024) and V — 1.45 (i.e., high inertia, low mass, and 
low pull velocity regime in the unsealed parameter space) 
as we decrease the dissipation parameter 7^ from 1.0 to 
0.001. We observe a wide variety of events, some of which 
are illustrated in the plots of the peel velocities (r) . For 



3 




25 

position 



FIG. 2: Plot of the peel velocity configuration for C/ = 
0.00788 and V" = 1.45. (a) A rugged peel front for 7„ = 0.1. 
(b) Illustrates a stuck-peeled configuration for 7u = 0.01. 



instance, all spatial points peel together for 7,1 — 1.0. As 
we decrease 7„ to 0.1 keeping other parameters fixed, 
the contact line profile becomes rugged even though all 
points peel nearly at the same time as seen in Fig. |2Ia). 

Intuitively, high 7„ implies that velocities of neighbor- 
ing points are coupled strongly and hence are not allowed 
to follow their local site dynamics. Thus, the total dissi- 
pation i?^g(T) = '^luY.ti.^i+i- ^if is vanishingly smaU 
when peeling is coherent. In contrast, for low 7„ [ say 0.1 
as in Fig. El^a)], the coupling between neighboring veloc- 
ities is weak and the local dynamics dominates. This 
means more ruggedness and hence higher dissipation. 

As 7i, is decreased to 0.01, the peel front exhibits two 
types of configurations depending on whether the sys- 
tem is on AB branch entirely, or partly on CD and AB 
branches of 0. When on CD branch, the ruggedness is 
substantially higher than that for 7^ — 0.1 (Fig. Ela)). 
Once the peeling process starts, the peel front breaks 
up into regions of stuck and peeled segments as shown 
in Fig. Elb). This configuration results from the orbit 
jumping between the low and high velocity branches of 
A typical phase plot of Xf versus vf is shown in 
Fig. Efa) for i = 25 with other points differing only in 
phase. Thus, as the phase difference along the peel front 
builds up to a value equal to the phase difference between 
the orbits that are in the stuck and the peeled states, the 
stuck state changes to a peeled state or vice versa. 

Now consider the dynamics for higher mass (C/ — 
0.788, vc = 0.024) and V = 1.45 as we decrease 7„ = 1.0 
to 0.001. For this entire range of 7u, the peel front dis- 
plays stuck-peeled segments for all times as shown in Fig. 
|3Jb) . There is a dynamic equilibrium between the peeled 
and stuck segments with the segments that are stuck at 
some instant getting unstuck at another instant and vice 
versa. Further, as is clear from Fig. |Sl[b), the average 
of the velocity jumps along the peel front is smaller than 
the low mass case (compare Fig. EJ. Concomitantly, the 
number of stuck segments increases with each stuck seg- 
ment having a only few stuck points better illustrated in 
an instantaneous plot oivf — l shown in Fig. Ol^c). More- 
over, from Fig. Ofc), it is clear that even the points that 
are in the stuck state are barely stuck. Further, it is clear 
that the orbit spends considerable time around the max- 
imum which is the critical peel value [Fig. Eld)]. Thus, 
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FIG. 3: (a) The phase plot for i = 25 for C/ = 0.00788 and 
V^" = 1.45, and 711 = 0.01. (b) Plot of 'an edge of peeling' con- 
figuration for Cf = 0.788 and V" = 1.45, and 7„ = 0.01. (c) 
The corresponding instantaneous plot of 'an edge of peeling' 
configuration. Dashed line represents the critical peel velocity 
'^rn.ax ~ 1- (d) The Corresponding phase plot for i = 25. Bold 
lines in (a) and (d) represent 4i{v"). 



FigsOfb) andini^c) correspond to a verge of peehng state. 
The 'edge of peeling' picture remains unaltered with time 
even though the stuck points themselves change. The 
largest Lyapunov exponent is 0.15 ( for N — 50) and 
hence this state is spatio-temporally chaotic [Fig. 01(a)]. 

In experiments, the nature of the AE signals changes 
from burst type to continuous type as the pull velocity 
is increased. In the model, the rate of dissipated energy 
i?*g = —dEae/dr represents the AE bursts. We have 
studied the statistics of as we increase the pull ve- 
locity keeping the tape mass low (C/ = 0.00788, t;c = 

0. 0024). As in experiments, for small 7„, at low veloci- 
ties, we find that R^^ exhibits bursts followed by a qui- 
escent state as shown in Fig. EJb) which is similar to 
Fig. 4(a) of Ref. 3] (for the AE amphtudes). [Fig. I^b) 
corresponds to Fig. I^b).] In contrast, for high pull ve- 
locities and low mass, and for a range of 7,j, R^^ exhibits 
continuous bursts as shown in Fig. ^c), which is again 
seen in experiments (Fig. 4(b) of Ref. ((3)- High mass 
and low pull velocity also exhibits continuous bursts. 

Denoting Efi to be the amplitude of i?^g (i.e., from a 
maximum to the next minimum), for high pull velocities 
and low tape mass ( Cj = 0.0078, Vc = 0.0024, = 
5.93), we find that the distribution of the magnitudes 
D{Er) shows a power law for all values of 7„ investigated, 

1. e., D{Er) - E]^"'^. Further, D{Er) shown in FigUd) 
exhibits two distinct scaling regimes as for case of the 
distribution of the AE amplitudes {A) in experiments. 
The value of toe 0.6 for the small amplitude regime 
while that for large amplitudes that has a substantial 
scatter is about 1.8. The corresponding exponent values 
are ~ 0.3 and 3.2 Using the fact that energy 
E oc A^, we get niE = (1-|-toa)/2. Inserting the values of 
toa, we get the corresponding exponents to be ttie 0.65 
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FIG. 4: (a) Spatio-temporal chaotic plot for C/ — 
0.788, V" = 1.45, and ju = 0.01 (color on line, low high 
peel velocity: red yellow green blue — + pink), (b) 
Plots of Ri^{r) vs. time r for C/ = 0.00788, V" = 1.45, 7„ = 
0.01. (c) for Cf = 0.00788, V' = 5.93, 7„ = 0.01 and 

(d) the corresponding distribution D{Er) of the amplitudes 
Er showing two scaling regimes. 

and tue ~ 2.1 which are close to the values predicted 
by the model considering the scatter for the latter. In 
contrast, for high mass and low velocity case, we find a 
single scaling regime with an exponent ~ 0.69. 

Thus, several qualitative features of the peel front dy- 
namics observed in experiments are reproduced by the 
model. For example, the characteristic features of the 
AE signals observed in experiments, namely, noisy AE 
bursts for low pull velocity changing over to continuous 
bursts at high pull velocity is reproduced. For high pull 
velocities ( low tape mass), D{Eb) exhibits two scaling 
regimes. However, comparison with experiments is made 
difficult due to the paucity of quantitative results except 
for the values of the exponents in the two scaling regimes 
which is in reasonable agreement with the model. Even 
so, our study suggest that if one wants a smooth peeling, 
one should peel at low velocity using high viscous dissi- 
pation. Significantly, the power law is seen at high pull 



speeds and thus unlike self organized criticality. 

The power law statistics for high pull velocities arises 
as a competition among the time scales due to inertia of 
the tape, dissipation and imposed velocity which is small 
at high leaving very little time for internal degrees of 
freedom to relax. The nature of peel front ranges from 
synchronous peeling for large 7„ to rugged type for small 
7tj. The 'stuck-peeled' configuration is qualitatively sim- 
ilar to the inhomogeneous peel fronts observed in exper- 
iments 1^ as also to the thin viscous film interface [l^. 
Interestingly, the 'verge of peeling picture' of the peel 
front (Fig. |31 b,c) is similar to the edge of unpinning 
picture of dislocations in the Portevin Le Chatelier effect 
[l4 | . This is one of the few cases where the the power law 
emerges purely from deterministic dynamics (see 

Here it is worth commenting on the assumption that 
the integrand of Eq. 1 vanish at each point y which is 
valid when L » b and when shear modulus fcg is small. 
In principle, one should have a long range term of the 
form kc jl^[u{y) — u(y')]'^dydy' /2\y — y'\. An equilibrium 
calculation with u{y) defined at one end shows that the 
shear strain energy is less than one percent compared to 
the total even for small b/L = 0.2. This lends support 
for the PE term used. We have also carried out numer- 
ical calculations by retaining this term. For small fcfj, 
results are not affected as it should be expected. How- 
ever, for relatively high values of ka, the solutions that 
were smooth breakup into stuck-peeled configurations. 

Finally, we note that properties of adhesive glue has 
been included in the model in an indirect way through 
the peel force function ( and low effective spring constant 
of the peel front due to adhesive glue) and that of AE 
through the Rayleigh dissipation function. We state that 
while the model recovers most dynamical features of peel- 
ing, issues that depend critically on the finite thickness 
of the adhesive material (for instance, fibril formation) 
cannot be addressed within the scope of the model. 
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